Bose-Einstein condensation in shallow traps 
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In this paper we study the properties of Bose-Einstein condensates in shallow traps. We discuss 
the case of a Gaussian potential, but many of our results apply also to the traps having a small 
quadratic anharmonicity. We show the errors introduced when a Gaussian potential is approximated 
with a parabolic potential, these errors can be quite large for realistic optical trap parameter values. 
We study the behavior of the condensate fraction as a function of trap depth and temperature and 
calculate the chemical potential of the condensate in a Gaussian trap. Finally we calculate the 
frequencies of the collective excitations in shallow spherically symmetric and ID traps. 

PACS numbers: 03.75.Fi, 32.80. Pj, 03.65.-w 



I. INTRODUCTION 

An atomic Bose-Einstein condensate is an excel- 
lent tool for studying quantum many body phenomena, 
such as collective excitations. The particle interactions 
in these condensates are weak and thus even quantitavc 
agreement between theory and experiments can be found. 
Recently a sodium condensate was trapped in an optical 
dipole trap Usually the potential of such trap is ap- 
proximated by the parabolic potential, but the shallow 
nature of the optical dipole trap makes it clear that for 
suitably large condensates, the anharmonicity of the po- 
tential has to be taken into account. 

In this paper we study how the shallowness of the po- 
tential changes the condensate properties. To enable an 
analytical approach we assume a Gaussian potential || , 
but we expect that the qualitative features are valid also 
for other types of shallow traps. Bose-Einstein conden- 
sation in power-law potentials has been studied previ- 
ously H , but usually the potential has been assumed to 
be parabolic or absent. In Sec. [n] we calculate how the 
condensate fraction behaves as a function of temperature 
and how the parabolic approximation may underestimate 
the BEC transition temperature considerably. In Sec. 
[II we calculate the chemical potential and estimate how 
many condensate particles can be trapped. Frequencies 
of the collective excitations are calculated in Sec. |v| and 
some concluding remarks are given in Sec. 

II. CONDENSATE FRACTION 

Particle interactions have a dramatic effect for parti- 
cle density distribution, chemical potential and collective 
excitations of the condensate, but the transition temper- 
ature and the condensate fraction can be accurately cal- 
culated with the ideal gas model . If there is a repulsive 
interaction between the trapped atoms, and the chemi- 
cal potential is almost the same as the trap depth, we 
expect changes to the ideal gas results, but for now, the 
simple ideal gas model will suffice. In a harmonic trap 
the transition temperature T c is given by 
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where io a is the trap frequency in direction a and N is the 
particle number. The condensate fraction is then given 

by [I 
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condensate fraction can behave in a qualitatively different 
way. Using a phase space density H 
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where H is the total energy, we can calculate the number 
of thermal particles as 
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Our trap has a finite depth and therefore we must intro- 
duce an appropriate r dependent cut-off for the kinetic 
energies. After a short calculation we see that the num- 
ber of thermal atoms is given by 
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where s max = exp(-x 2 ). 

When the trap has only a few eigenstates, the con- 
tinuum approach fails and we should model the system 
using a discrete spectrum. If the potential is spher- 
ically symmetric and is approximated as a parabola, 
then the number of different energy levels is roughly 
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Vi/hoj = yVima 2 /h 2 . This provides a lower limit to 
the number of energy levels, since for a Gaussian trap 
the separation between adjacent levels becomes smaller 
as the energy increases (approaching zero as the energy 
approaches the trap depth). If the trap depth is very 
small, Vi Rt 9 nK and a = 15 /jm the trap has roughly 
ten energy levels. We expect that for deeper traps the 
continuum approach should give accurate results. 

Let us now consider consider 10 6 sodium atoms in a 
trap with {(Jx&y&z) 1 ^ = 15/im. As the first example 
we fix the temperature to the value T = 300 nK and 
vary the trap depth. In Fig. |l| we show the resulting 
condensate fraction and compare it to the result we get 
by approximating the potential as parabolic. At large 
trap depth the condensate fraction is well predicted by 
the parabolic model, but at smaller trap depths the be- 
havior changes qualitatively. Parabolic model predicts 
that the condensate fraction should vanish, together with 
the trap frequency. This is obviously not the case. As 
the trap depth becomes smaller, the condensate fraction 
takes a minimum value after which it approaches unity. 
The minimum condensate fraction is achieved when the 
trap depth is about the same as the temperature. This 
happens because in a shallow potential there is simply no 
room for thermal atoms. At the ultimate limit the po- 
tential would have only one bound state and there would 
not be any states accessible to thermal atoms. We also 
see that at larger trap depths the condensate fraction is 
smaller than the result for the parabolic trap; a sensi- 
ble result since anharmonicity makes the trap more open 
than a purely parabolic result, thus reducing the effective 
trap frequency and lowering the critical temperature. 

As a second example we fix the trap depth V\ = 1 /J.K 
and vary the temperature. In Fig. ^| we show the con- 
densate fraction and compare it to the parabolic result. 
Again it can be seen that the behavior is dramatically 
different in a shallow potential. The condensate fraction 
is considerably larger over a wide range of temperatures 
and the critical temperature is about four times larger 
than the value predicted by the parabolic model. At 
very small temperatures (T < 240 nK) the openness of 
the Gaussian trap is reflected in condensate fractions, 
which arc smaller than the parabolic trap predictions. 
But this effect is so small that it is not visible in Fig. ||. 



III. CHEMICAL POTENTIAL VERSUS 
PARTICLE NUMBER 

Properties of pure condensates are well described by 
the Gross-Pitaevskii (GP) equation 

h 2 

v 2 * + y* + Nu \y\ 2 y = /i*. (7) 

2m 

Here <3/ is the condensate wavefunction, m is the atomic 
mass, V is the trap potential, N is the number of atoms 
and Uq = Airh 2 a s /m, where a s is the s-wave scattering 



length. When the particle number N is large, the ki- 
netic energy becomes small compared to trapping and 
the atomic interaction energies. In the Thomas-Fermi 
(TF) approximation the kinetic energy is ignored and we 
get an analytical result for the condensate wavefunction 



\V(x,y,z)\ 2 = 



1 



NU 



(p-v), 



(8) 



when the R.H.S. is positive and zero elsewhere. The 
wavefunction is normalized to unity, so by integrating 
Eq. (||) we get a formula that relates the number of par- 
ticles to the trap geometry and to the chemical potential 
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where error function is defined as 
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and x = fi/Vx- It is clear that the result converges only 
if < fi < V\. When \x = V\ we cannot add any more 
particles to the condensate, as these extra particles can 
not be trapped. From this condition we get the maximum 
condensate particle number in a Gaussian trap as 
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It is instructive to calculate this number for a sodium 
(a s = 2.75 nm JhJ) condensate with the reasonable 
trap parameters Vi = 4/iK, a x = cr y = 3 /im, and 
a z = 38 a x Q. The maximum number of condensate 
atoms is then very large, 
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For this type of trap the maximum condensate number 
density would be about 



Vx/Uo « 5.5 



10 15 cm" 3 , 



(13) 



and the three-body decay would limit the condensate life- 
time considerably. It should be noted that for conden- 
sates with the maximum number of atoms, the maximum 
density depends only on the trap depth and interaction 
parameter and is independent of the trap geometry. This 
is a general result and applies to any shallow trap as long 
as the TF approximation is valid. 

Equation (J9J) is somewhat awkward and it is useful 
to derive a simple approximation for it. The chemical 
potential is often much lower than the trap depth and 
we can then expand Eq. (^) around the small parameter 
x = fM/Vi. Keeping terms up to order x 5 / 2 we get the 
result 
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which corresponds to approximating the potential as 
parabolic and reduces to the familiar formula for 
parabolic potentials when we notice that the trap fre- 
quencies are related to a and Vi by ui = y/Vi/ma 2 . 
Keeping terms up to order x' 7 1 2 we get an equation 
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The solution to this equation is expected to be close to 
the first approximation so we can set x = xq + 5x (xq = 
Ho/Vi) and solve for 8x. The chemical potential is then 
given by 
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We have compared this formula to the exact numerical 
result for a condensate in a spherically symmetric trap 
and noticed that it is an excellent approximation even for 
relatively large values of x. For example when x 0.5 
Eq. (jll]) is exact with a relative accuracy better than 
1CT 3 . 

Let us calculate how large this shift in /i is for some 
reasonable parameters. The condition xq <C 1 implies 
that 
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must be satisfied. For a sodium condensate this means 
that 



W^ >>2 , 10 -51 Jm 3 i 



(18) 



a condition that is not difficult to achieve. Choosing 
reasonable parameters a x — a y — 3 fim, a z — 38 a x , 
Vx = 4.0 fiK and N — 10 6 || we see that Eq. @ predicts 
the shift in the chemical potential to be about —6 %. This 
shift is quite large, and it might be necessary to take this 
shift into account, when studying condensates in realistic 
optical dipole traps. 



IV. COLLECTIVE EXCITATIONS 

In the following subsections we calculate the frequen- 
cies of the collective excitations in spherically symmetric 
and ID traps. We aim at simple analytical results that 
clarify the role of trap anharmonicity and therefore we do 
not study anisotropic traps. The collective excitation fre- 
quencies for anisotropic traps can be solved numerically, 
but analytic results are exceedingly difficult to obtain. 



A. Spherical trap 

Collective excitations for a spherically symmetric 
parabolic trap with trap frequency f2 have been calcu- 
lated in the TF limit analytically pi . These excitations 
have the form 5n(r) = P 2nr (r / R)r r Yi m (9, (f) where P 2nr 
are polynomials of degree 2n and the dispersion law is 
given by the formula 
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In an anharmonic trap these frequencies will be shifted, 
but it is not known by how much. We assume spherically 
symmetric potential 
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where AV(r) — ar A . In particular, if the exact trapping 
potential has the Gaussian shape 

V exact (r) =Vx(l- eM-r 2 /2a 2 )) , (21) 

we can approximate it with a Taylor series and get 
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We choose the unit of length to be L = ^jTi/mVl and 
unit of time as r = I/O. In these new units the CP- 
equation becomes dimensionless 
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with the dimensionless interaction parameter a — 
ATra s /L and potential 



V(r) = ^r 2 - er 4 , 



(25) 



where e = hil /8Vi. 

Following Stringari [ll| ] we write the wavefunction in 
terms of phase and modulus, 



m(r,t) = ^n{f,t)/Ne i * { - f < t \ (26) 

where n is the density and velocity is fixed by the relation 

v(f,t) = (h/m)V4>. (27) 

The GP equation is equivalent with two coupled equa- 
tions 
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In the Thomas-Fermi limit the kinetic pressure term can 
be neglected in the equation for the velocity field, which 
then becomes 
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If we set d a == o {ot — x, y, z) and the chemical potential 
is sufficiently small, the stationary solution of Eqs. ( |2"g| ) 
and ( |30"| ) coincide with the Thomas-Fermi wavefunction 
studied in the previous section. Linearizing Eqs. (^8|) and 
( |30| ) by setting n = Uq + 6n(r)e tuJt gives us an equation 
for the time-dependent solutions 

V • [c 2 (r)V6n] + u 2 8n = 0. (31) 

Here uq is the stationary solution and mc 2 (r) = fx — 
V(r). We now calculate the corrections to the collective 
excitation frequencies as we move from a parabolic trap 
into an asymmetric trap. We proceed as in the Ref. Ill . 

In our case mc 2 = (i - ln §-r 2 + -^r 4 . The edge R of 
the condensate is defined by 
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With this definition the equation 
Vi 



R 2 



R 4 



(32) 

l|) takes the form 

Y,)n ). (33) 



+2 {n) 6n = °- 

We choose the dimensionless distance as x — r/R and 
frequency u) 2 D = 2(ui/tt) 2 . If we postulate a solution 5n = 
P(r)Yim(6, 4>) and define a dimensionless parameter f3 = 
■^s = [i/2Vi+0((/i/Vi) 2 ), we obtain an equation for the 
radial part 
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We can consider the terms proportional to (3 as a per- 
turbation and try to find a first order solution to the 
equation with form 



HP + AHP 
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and P is the known solution without anharmonicity. The 
frequency ui of the collective excitation is then given by 



where 
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and Uq = 2n 2 



2n r l + 3n r + I are the unperturbed fre- 
quencies. As the result is perturbative and expected to be 
relevant only for the lowest excitations we will give results 
only for surface excitations (n r = 0) and the breathing 
mode (n r = 1 and I — 0). 

For surface excitations (P = x l ) we get the result 



21 + 5 



and for the breathing mode (P — 1 — | x 2 ) we get 
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In Fig. |3| we com par e the exact numerical solution, to 
the approximation ( |39| ) , and to the Thomas- Fermi result 
for the modes (n r = 0, 1 = 1) and (n r = 0,1 = 2). We 
assume 10 6 sodium atoms in a trap with a = 15 /im and 
vary the trap depth. We see that our analytical results 
match the exact numerical values quite accurately. Since 
the TF approximation was always justified, the numeri- 
cal values for the parabolic case are well predicted with 
the well known analytical results. If, on the other hand, 
number of particles would have been less, we would ex- 
pect (n r = 0, 1 = 2) mode frequency to decrease signifi- 
cantly with increasing trap depth, approaching the value 
y/2 asymptotically. This trend is opposite to the results 
from a fully Gaussian trap and is a clear indication of the 
failure of the parabolic approximation. 

In Fig. |] we show the corresponding comparison for the 
lowest compressional mode, the breathing mode. Again 
it can be seen that Eq. (^) gives an accurate approx- 
imation to the exact result and that corrections to the 
parabolic model are noticeable. 



B. ID trap 

Spherical symmetry, discussed in the previous subsec- 
tion, is a special case not necessarily easily obtained in 
experiments. To enable analytical approach and to get 
a feeling about possible changes due to different trap ge- 
ometries we will now calculate the collective frequencies 
also for anharmonic one dimensional trap. The potential 
is given by 



V(x) = 



(41) 
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From the TF limit we will get an equation relating chem- 
ical potential and condensate boundary 
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Using this result and continuing in the same manner as 
in previous subsection we get an ID analogue of Eq. (p4) 
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where (3 = eR 2 /mu 2 . If the anharmonicity results from 
a Gaussian potential, then (3 = R 2 /8a 2 . When (3 — 

the solutions have the form Sn = x P Y2k=o akX > where 
p = or 1. If p = solutions are even and excitation 
frequencies are given by ux^ = 1/2 • N(N + 1). If p = 1 
solutions are odd and corresponding frequencies are given 
by uj 2 n = 1/2 • (N + 1)(N + 2). For both even and odd 
solutions N must be an even integer. 

For the lowest excitations we calculate the shifted fre- 
quencies uo 2 (N = 0,p = 0) = 0, u 2 (N = 0,p = 1) = 
1 - 4(3, uo 2 (N = 2,p = 0) = 3 - ^(3 and u 2 (N = 2,p = 
1) = 6 — ^-(3. Especially we see that the corrections are 
always proportional to the dimensionless parameter (3, in 
the same way as in the spherically symmetric case. It is 
to be expected that the calculated shifts will give correct 
order of magnitude estimates even for more complicated 
trap geometries. 



V. CONCLUSIONS 

In this paper we have demonstrated several ways by 
which the shallow anharmonic potential changes the 
condensate properties. We have limited our discussion 
mainly to Gaussian potentials, but most of the results 
can also be applied to traps having a quadratic anhar- 
monic term. Naturally, as an order of magnitude esti- 
mates, these result should be applicable also for other 
types of shallow traps. It seems that in optical dipole 
traps the errors introduced by approximating the poten- 
tial as harmonic can be quite large. Harmonic model 
can predict the condensate fraction poorly and with rea- 
sonable parameters the corrections to the chemical po- 
tential and collective excitation frequencies can be sev- 
eral percent. In an optical dipole trap the spin degree 
of freedom is not necessarily frozen and the condensate 
should be described with a multicomponent spinor wave- 
function |l2| Q . Nevertheless, the results in this paper 
might prove to be useful also in studies of spinor conden- 
sates. 
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FIG. 1. The condensate fraction as a function of trap 
depth for 10 6 sodium atoms, when (cr a; (T y cr 2 ) 1 ' /3 = 15 fim and 
temperature T — 300 nK. The solid line is the numerically 
calculated value using Eq. (^|) and the dashed line is the result 
we get by approximating the Gaussian potential as parabolic. 



FIG. 3. Frequencies of the collective excitations with 
(n r = 0, I = 1) and (n r = 0, 1 = 2) as a function of trap 
depth, for N = 10 6 sodium atoms and a = 15 /im. The solid 
line is the exact numerically calculated value, the dashed line 
is the numerically calculated result for the parabolic approx- 
imation, and the dot-dashed line is based on Eq. The 
exact result and the approximation (^9|) give almost the same 
results for most values of Vi. 




FIG. 2. The condensate fraction as a function of temper- 
ature for 10 6 sodium atoms, when (cr^cr^crz) 1 ^ 3 = 15 /im and 
trap depth Vi = 1 pK. The solid line is the numerically cal- 
culated value using Eq. (^) and the dashed line is the result 
we get by approximating the Gaussian potential as parabolic. 
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FIG. 4. Frequency of the breathing mode (n r = 1,1 = 0) 
as a function of trap depth for N = 10 6 sodium atoms and 
a = 15 fim. The solid line is the exact numerically calculated 
value, the dashed line is the numerically calculated result for 
the parabolic approximation, and the dot-dashed line is based 
on Eq. @. 
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